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Abstract 

We analyze the existed solar neutrino experiment data and show the allowed regions. The 
result from SNO's salt phase itself restricts quite a lot the allowed region's area. Reactor neu- 
trinos play an important role in determining oscillation parameters. KamLAND gives decisive 
conclusion on the solution to the solar neutrino puzzle, in particular, the spectral distortion 
in the 766.3 Ty KamLAND data gives another new improvement in the constraint of solar 
MSW-LMA solutions. We confirm that at 99.73% C.L. the high-LMA solution is excluded. 
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I Introduction 

The electron neutrinos emitted from the sun disappear somewhere when they travel to the earth. 
This is the famous solar neutrino deficit, which is the almost forty years' "Solar Neutrino Problem" . 
There were many attempts to solve this puzzle during the years. Some of them were tried to 
modify the solar model in order to give a lower original neutrino flux, which conflict the energy 
spectrum provided by the 4 first-generation experiments: Homestake, Sage, Gallex and Kamiokande 
El H|- Recent experiments have shown that the solar neutrino oscillate by v e — ► inside 
the sun via MSW conversions. This was proven by the Sudbury Neutrino Observatory (SNO) 5 , 
and it was confirmed by the reactor experiment KamLAND 0. The former experiment detects v e , 
Ve J r l/ fjL + v T and v e + 15%(^„ + u T ) three quantities on earth, which correspond to CC, NC and 
ES interactions respectively; KamLAND observes v e — > u x neutrino oscillation channel. 

II Solar neutrinos 

The solar neutrino puzzle was solved by the neutrino oscillations u e — » z/„ r inside the sun via MSW 
conversions. This was proved by the Sudbury Neutrino Observatory (SNO) in Canada. And it was 
confirmed by the laboratory base line experiment KamLAND in Japan. 

SNO is a 1000 ton heavy water Cerenkov detector mainly measuring S B solar neutrinos. It 
consists of nearly 9450 photon-multiplier tubes and light concentrator units arrayed on a geodesic 
support structure, with light water surrounding the spherical acrylic vessel containing the Z^O.The 
first phase of SNO data is from the pure D2O. After that the experimenters add up NaCl (salt) 
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to enhance the NC events rates. This is called the second phase or "salt phase". 
In analysis of the solar oscillation data [7j, we use the x 2 defined as: 



X® — Xlgen + XSK + XSNO 



(2.1) 



where Xigen s ^ an< ^ s f° r Chlorine and Gallium experiments. To calculate each individual chi square 
in the right hand side of eq. (|2.1|) . we use the so called covariance approach: 

N 

Xcovar = (Rn P ~ -R-n )[°"nm] {R"mF ~ ^"m )■ (^-2) 

n,m=l 

Here i?n Xp and correspond to experimental result and theoretical value for the n-th data 

point. N=2, 34,44 are for x\ ge n > Xski Xsno respectively. For getting a i?* °, the important step 
is to calculate the v e survival probability. We have used three methods to check its consistency: 
the Parke formula jHj; the modified semi-analytic formula in 0; and the completely numerical 
propagation. We found that the second way is the best, considering both the calculation precision 
and the computer CPU hour. 

The covariant matrix of squared error cr^ m can be written as 



/i n is the uncorrelated error of the n-th detected quantity for both the experiment and the theory 
(such as the statistic uncertainty, or other uncertainties which affect only one detectable value), 
and is the correlated systematic error caused by the k-th correlated original error (the original 
error may be the spectrum uncertainty, or the energy resolution uncertainty, etc.). For a detectable 



K 




(2.3) 



k=i 
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value i?n Xp , we can say that its uncertainty (R^ p — R^ ) is in the range of ±/i n ± ± ± ... ± c^. 

For Xigem ^ ne correlated errors are the 12 SSM uncertainties (the cross section factors 
Sn, "Si3; S'34, £1^4, S17; the 7 Be capture cross section Ce e ; the solar luminosity; metallicity Z/X; 
age; opacity; the element diffusion and Sh ep ), whose relative uncertainties AlnXj. determine the 
correlated uncertainties of the neutrino fluxes &f SM through the logarithmic derivatives 

With this uncertainties and the matrix given in (|2.3|) (|2.4j) . we can calculate the fractional 
uncertainties of the SSM neutrino fluxes: 



A In $f 5M 
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\ k=l 



For Xsno an( i Xsio the om y effective SSM uncertainty is the B neutrino flux uncertainty since 
the hep neutrino flux is too small. So, we uses free S B flux method ^01; i-e., we float the flux near 
the central value given in paper within la uncertainty range, then find out which flux gives 
the minimal \ 2 ■ The \ 2 formula changes a little as 



N h - f V 2 



n,m=l 

Thus to describe a\ m in these two experiments, we don't need to consider the SSM uncertain- 
ties. The remaining original uncertainties that can affect x\k are ^ ne spectrum shape error, 
the energy scale uncertainty, the resolution uncertainty and an overall SK systematic offset un- 
certainty. For SNO, the remaining original uncertainties are S B spectrum shape error, the energy 




Figure 1: The 90%, 95%, 99%, and 99.73% C.L. allowed regions in the Amf 2 - tan 2 6 12 plane after 
SNO experiments. Left corresponds to the result of SNO's first phase, the right is the result of 
SNO's salt phase case. 

scale uncertainty, the resolution uncertainty, the vertex reconstruction uncertainty, the cross section 
uncertainty, the neutron capture uncertainty, the neutron background uncertainty, and low energy 
background uncertainty [S]. 

Using this method we scan the solar neutrino tan 2 812 — Amj 2 parameter space, neglecting the 
tiny correction from #13 and Am 2 3 . The allowed regions are shown in figs. ^ 

III Manufactured neutrinos 

At present the only detected laboratory neutrinos which are sensitive to Am 2 2 — tan #12 parameters 
are reactor neutrinos. KamLAND is a first experiment to deal with these solar neutrino parameters 
completely from manufactured source instead of from the sun. This reactor neutrino experiment 
convinces scientists that the Large Mixing MSW is the solution of the solar neutrino problem 

[running. 
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In a reactor, anti neutrinos are released by radioactive isotope fission; the total neutrino spec- 
trum is a rather well understood function of the thermal power W, the amount of thermal power 
Wi emitted during the fission of a given nucleus, and the isotopic composition of the reactor fuel 

The index i of /, stands for 4 isotopes such are 235 U, 238 U, 239 Pu and 2Al Pu. The (dN/dE) is the 
energy spectrum of the fissionable isotope, it can be parameterized by the following expression |14| 
when E v > 2MeV: 

e a +a 1 E v +a 2 El /g 2) 



dE v 

the coefficients ai depend on the nature of the fissionable isotope. KamLAND is a scintillator detec- 
tor, where electronic anti neutrinos are detected by free protons via inverse j3— decay reaction |14j. 

V e + p — > e + + n (3-3) 

in the limit of infinite nucleon mass, the cross section of this reaction is given by cr(E u ) = kE e +P e + , 
where E e +, P e + are the positron energy and momentum respectively and k can be taken as k = 
9.55 x 10 -44 cm 2 /MeV 2 . The anti neutrino events are characterized by the positron annihilation 
signal and the delayed neutron capture sign [Sj. 

From the reactor to the detector, massive neutrinos oscillate on the way and change their flavor 
composition to a certain extent. The anti neutrino, V e , can oscillate to other flavors via Amf 2 and 
Am 2 n. For KamLAND experiment, the distance is in the range of hundred kilometers. Due to the 
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Figure 2: The combined LMA region of solar and KamLAND data before conference Nu2004. The 
difference between left and right shows improvement from SNO's salt phase result. The confidence 
levels are the same as what in fig. ^ 

tiny value of #13 and the value of Amf 3 , contribution to oscillation from the biggest neutrino mass 
scale 7713 gives a small averaged effect thus we can reduce it to a two-flavor neutrino analysis: 

P(E„, L, 9 12 , Am? 2 ) = 1 - sin 2 (2tf 12 ) sin* (^^^) (3-4) 

The combined results of solar and reactor neutrino experimental data before "Neutrino 2004" 
are shown in figs. El 

Using the newest 766.3 Ton- Year of data published by KamLAND collaboration at conference 
Neutrino 2004, we find the tan 2 612 — Am^ allowed region with the energy spectrum in 13 bins. 
The result is shown in fig. El-left. The combined chi square is then x 2 = X® + Xkl with result 
shown in fig. -right; which confirms the Large Mixing Angle MSW solution to the solar neutrino 
problem. The best fit point we get is 

tan 2 9 = 0.4; Am^ = 8.3 x W^eV 2 , 
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Figure 3: Left is the new KamLAND allowed regions for 766.3 Ty data, together with the LMA 
solution contours (not filled in) by solar neutrinos themselves. Right is the combined LMA region 
of solar and KamLAND data. The confidence levels are the same as what in fig. ^ 

which is consistent with the previous papers ^Sji The day-night effect ^H] f° r solar neutrinos is 
then expected to be about two to four percent. 

Many other possible candidates for the solution of solar neutrino problem are excluded; sterile 
neutrinos are no more supported in this long-lasting puzzle. However, the most reliable candidate, 
the LMA MSW solution is becoming precise and more confident in terms of a much larger confidence 
level and much smaller allowed region. Thus we conclude that the solar neutrino problem is solved 
by large mixing MSW solution, with parameters shown in fig. [2 

In conclusion, using the new data of KamLAND experiment we conclude that the solar neutrino 
solution is in the large mixing MSW adiabatic allowed region, corresponding to (tan 2 6, Am 2 ) ~ 
(0.27 - 0.62, 7.2 x 10" 5 - 9.5 x 10~ 5 eF 2 ), at 4a level. Which is an amazing small region in the 
parameter space with such a high confidence level. Our result with best fit point at (tan 2 9, Am 2 ) ~ 
(0.4, 8.3 x 10~ 5 eV 2 ) is in good agreement with previous studies. Predictions of the MSW adiabatic 
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conversion solution of the neutrino energy spectrum, the day-night effect, the seasonal variation, 
the ratio of CC and NC current events are all consistent with the solar data. At present, none signal 
indicates other solutions rather than the MSW adiabatic solution, together with the standard solar 
model. 
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